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Abstract 

The electrocaloric effect (ECE) in normal and relaxor ferroelectrics is investigated in the frame- 
work of a thermodynamic approach based on the Maxwell relation and a Landau-type free energy 
model. The static dielectric response of relaxors is described by the spherical random bond-random 
field model, yielding the first Landau coefficient a = a(T), which differs from the usual expression 
for ferroelectrics. The fourth-order coefficient b is treated as a phenomenological parameter, which 
is either positive or negative due to the anisotropy of the stress-mediated coupling between the 
polar nanoregions. When b < 0, the maximum ECE in a relaxor is predicted near the critical point 
in the temperature-field phase diagram, whereas in a ferroelectric it occurs at the first order phase 
transition. The theoretical upper bound on the ECE temperature change is estimated from the 
values of saturated polarization, effective Curie constant, and specific heat of the material. 

PACS numbers: 77.70. +a,77.80.Jk,77.84.-s 
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I. INTRODUCTION 



Several articles have recently focused on the electrocaloric effect (ECE) in ferroelectrics 
and related materials,- - - which bears analogy with the well known magnetocaloric effect 
(MCE).— £ Here we investigate the mechanisms of ECE in relaxor ferroelectrics, to be referred 
to as relaxors and normal ferroelectrics (or ferroelectrics), and discuss the specific features 
of these two groups of materials. In particular, we will discuss the possibility of achieving 
a giant ECE in bulk inorganic relaxors and ferroelectric materials as well as in organic 
polymers. These systems offer the prospect of practical applications, such as miniaturized 
and energy efficient cooling devices, without the need for large electric currents commonly 
associated with the MCE. 

A crucial physical quantity in ECE is the change of entropy of a polar material under 
the application and/or removal of an external electric field. For example, when the electric 
field is turned on isothermally, the elementary dipolar entities in the system will become 
partially ordered and the entropy will be lowered. The entropy lowering of the dipolar 
subsystem is then compensated by an increase of the temperature of the total system, 
which characterizes the ECE. The degree of lowering depends on the number of statistically 
significant configurations in the initial and final states of the system, as well as on the size 
of the average dipole moment and the volume density of dipolar entities. Other factors may 
also play a role: If the system undergoes a first order phase transition under the action of 
external electric field, the entropy will be enhanced on crossing the borderline between the 
two phases, resulting in a larger ECE. The line of first order transition points terminates at 
a critical point where the transition becomes continuous,- and it will be of special interest 
to investigate the behavior of ECE in the vicinity of the critical point. 

Estimates of the ECE can be made on the basis of thermodynamic Maxwell relations 
using the measured heat capacity and the field and temperature dependence of the dielec- 
tric polarization. From the theoretical point of view, a central problem is how to make 
predictions about the temperature and field dependence of ECE. As a first step, one needs 
to develop an appropriate phenomenological and/or mesoscopic model, which incorporates 
the specific physical features of the systems. Here we will make use of the standard Landau 
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phenomenological model, which can be applied to both relaxors and ferroelectrics with the 
corresponding choice of Landau coefficients. These in turn can be derived from the meso- 
scopic model of the material under study. In the case of relaxors, the mesoscopic model of 
choice is the spherical random bond-random field (SRBRF) model, which is based on the 
concept of reorientable polar nanoregions (PNRs). 8 Thus we should be able to compare the 
ECE in relaxors and ferroelectrics, and determine the parameters, which control the ECE in 
these systems. Finally, using general principles of statistical thermodynamics we will discuss 
the existence of a theoretical upper bound on the ECE and argue that it satisfies a universal 
relation, which is, in principle, also applicable to MCE. 

II. THERMODYNAMIC APPROACH 

The temperature change of a polar system under adiabatic electric field variation from 
the initial value Ei = to final value Ef = E can be written in the form 1 



which follows from the thermodynamic Maxwell relation (dS/dE)T = (dP/dT)E involving 
the entropy density S(E,T) and the physical dielectric polarization P (in units of C/m 2 ). 
The volume specific heat at constant field is given by Ce = PCe- 

In deriving Eq. ([Q), one tacitly assumes that the fluctuations of polarization P(r) can be 
ignored and that P represents a thermodynamic variable given by the macroscopic average 
of P(f). Furthermore, it is implied that the system is ergodic, i.e., its response time much 
shorter than the experimental time scale. 

If the field and temperature dependence of P(E, T) is known from experiments, the 
integral in Eq. ([1]) can be evaluated, yielding an estimate for AT.-^ 

In model calculations, it seems convenient to change the integration variable in Eq. ([TJ 
from dE to dP(E). This is readily done by applying the thermodynamic identity 9 




(1) 




(2) 



with the result 




(3) 
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This expression is fully equivalent to Eq. ((T]), with the new integration limits given by 
P Q = P(0,T) and P = P{E,T). 

The partial derivative (dE/dT)p can be obtained from the free energy density functional 
F(P,T). Ignoring fluctuations of the order parameter P, we write F as a power series 

F = F + -aP 2 + -6P 4 + -cP 6 + EP. (A) 

2 4 6 V ; 

This has the standard form of a mean field free energy expansion with temperature dependent 
coefficients a,b,c, etc. 

Applying the equilibrium condition (dF/dP)x = 0, we obtain the equation of state 

E = aP + bP 3 + cP 5 + ---, (5) 

and the temperature derivative in Eq. (J3j) becomes 

(^) p = aiP + 6ljP3 + ClP5 + '"' (6) 

where a± = da/dT, b\ = db/dT etc. It should noted be that P = P(E,T) in Eq. is that 
solution of Eq. (JSJ), which simultaneously minimizes the free energy (T4|). 
The integration in Eq. fl3]) can now be carried out, yielding 



T 

AT= — 



(7) 



In passing, we note that Ce, in general, depends on the temperature; however, in writing 
down Eqs. flTJ and (j3J) the temperature dependence of the heat capacity had already been 
ignored. 

The expression in brackets is related to the change of the entropy density S = 
— (dF/dT)E- Using Eq. (jlj) we can write 

S = 5 + 5 1 (P). (8) 

The first term So = —dFo/dT is the entropy at P = 0. It contains the configuration entropy 
of dipolar entities, which depends on the number of equilibrium orientations Q, say, Q = 8 
for the < 111 > equilibrium case.— Thus we may expect that Sq ~ (N/V)khx(Q), N being 
the total number of dipolar entities such as PNRs in relaxors. The second term is given by 

Sl (P) = " (I ai ^ + \ h ^ + \ 01 ^ + " ' ) • (9) 



Back in Eq. (J7|), Sq cancels out and the ECE temperature change can be rewritten in the 
familiar forrn^ 

AT = ~AS P , (10) 

with AS P = Sx(P) - Si(P ). 

It should be noted that the values of all temperature-dependent quantities Pq, P, a\, b±, 
etc., on the r.h.s. of Eq. (j7]) are taken at the initial temperature T = T, and Cp at the final 
field value E. 

The coefficients a, b, c, ... can be expressed in terms of linear and nonlinear susceptibilities 
by formally inverting the relation (J5]) and writing P(E) as a power series in EM- In Landau 
theory, close to a second order phase transition one sets a(T) oc T — To, while b, c, ... are 
constants. Thus, a\ = const., and b\ — c\ — • ■ ■ — 0. This leaves only one nonzero term 
of the order 0(P 2 ) in Eq. (J7J). On the other hand, Xi oc |T — T | _1 and the nonlinear 
susceptibilities are also found to diverge when T — > T . Thus a formal inversion P(E) of 
Eq. ((SD in powers of Tv would lead to a poorly converging series. 

In the following we will apply Eq. (J7J) in order to discuss the predictions of Landau theory 
in two characteristic cases, namely, normal ferroelectrics and relaxors. 



III. SRBRF MODEL OF ECE IN RELAXORS AND FERROELECTRICS 

As already mentioned, in Landau theory of phase transitions in ferroelectrics, the coeffi- 
cients b,c, ... in Eq. (j4]) are assumed temperature independent and a oc T — T , where T is 
the Curie- Weiss temperature. When b > and E = 0, a second order transition occurs at 
T c = Tq. For b < and c > 0, a first order transition appears at a temperature T 1; given by 
the relation a(T) = 3b 2 /16c. Writing a = (T — T o )/(eo0), where 9 is the Curie constant, 
we find: 

T 1 = To + £ e4-- (11) 

lb c 

For E ^ 0, a critical point will be located at T^p, Tv^p, where^ 

^ 9 6 2 ^ 66 2 /3j&T , x 

T CP = T + e e , E CP = — \/ — . 12 

° 20 c 25c V 10c V ; 

Turning next to relaxors, we assume that the relevant elementary dipolar entities at 
temperatures around the dielectric maximum are polar nanoregions or PNRs. According 
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to the SRBRF model,- these PNRs are coupled through Gaussian random interactions Jjj 
("random bonds") and are subject to Gaussian random fields hi. In zero applied field, 
spontaneous long range order is suppressed (Pq = 0). This means that for a relaxor we 
can still use the free energy 01]), however, the coefficient a(T) must remain positive at all 
temperatures. Thus, for b > and E = there can be no second order phase transition. 
The explicit form of a(T) follows from Eq. fl5]), namely, 



dE 
dP 



= (e Xi)~\ (13) 

p=o 

where xi — [9P '/ '9E)e=o represents the (quasi)static linear field-cooled dielectric suscepti- 
bility. This can be derived from the SRBRF model of relaxors: 8 

= 9(1 ~ q) (14) 

Xl T-T (i- q y 1 1 

with 

VKEq 

The effective Curie constant is given explicitly in terms of the average squared dipole 
moments Qi of PNRs, namely, g 2 = Y2iLi 9i/N, and the average volume v = V/N associated 
with a PNR. Later, we will also introduce the saturation polarization P max = g/v, where 
g = Yji=i9i/N- For simplicity, we will henceforth neglect the difference between g and g, 
and write P max = g/v. 

The parameter can be determined experimentally from the asymptotic high temper- 
ature behavior of Xi-~^ The parameter To = Jo/k is defined via the average over the 
infinitely ranged random interaction [Jij] a „ = Jq/N. The spherical glass order parameter, 
< q < 1, is a measure of the degree of disorder. For P = it is determined by the cubic 
equation 

{kTf q ={J 2 q + A){q-lf. (16) 

Here, J 2 is proportional to the variance of the random bond distribution according to [Jfj] a v — 
(\Jij]av) 2 = J 2 /N , while A measures the correlations of quenched random fields, i.e., [hj\ av = 
and [hihj] av = 5ij A. If J < y/ J 2 + A, long range order will be suppressed (Pq = 0). Thus 
the relaxor state is characterized by three physical parameters: Jo, J, and A. Typically, 
in relaxors one finds A < J 2 , implying that random bonds are effectively much stronger 
than random fields. This allows PNRs to reorient collectively under the action of external 
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fields and relax towards equilibrium. In the opposite case, A J 2 , PNRs are would be 
trapped in a frozen static configuration of random fields and no characteristic low-frequency 
relaxor response due to PNR flipping could be observed. Here we will consider only the case 
A < J 2 . 

From Eqs. (I14p and ( II 3p we derive the coefficient a(T): 

a (T) = —7c { t~~~ ~ Tq \ ■ (17) 



e Q \l - q 

The linear susceptibility \i has been fitted to experimental data for the static field-cooled 
response in a variety of relaxor systems, from which the parameters of the model have been 
obtained.— r — 

Formally, the SRBRF model also yields explicit expressions for nonlinear susceptibilities, 
from which the coefficients h and c in the free energy (j3J) can be determined.— However, it 
has been shown earlier— that realistic values of these coefficients can only be obtained if 
the coupling between PNRs and lattice strain fluctuations is included. Several mechanisms 
for such a coupling have been investigated both at a mesoscopio^"^ 1 ^ and phenomenological 
level.— It has also been shown^^ that in real three dimensional systems strain coupling 
gives rise to anisotropy of the anharmonic terms in the free energy— Specifically, strain 
coupling may change the sign of the coefficient b and hence of the corresponding nonlinear 
susceptibility X3 f° r a given direction of the applied field. In the following, we will simply 
consider b and c as free parameters and discuss separately the cases b > and b < 0, while 
keeping c > 0. 

For cubic systems, Eq. (jlj) can be rewritten in a general form 20 

F = a x Y,Ff + a n J2 P i + a ^J2 P i P I + «m 

i+i i (lg) 

+ «i 22 P ^ P ? + P ") + <*™ P MPt ~ E E i p » 

where «i, a^-, are so-called dielectric stiffness coefficients with i,j,k = 1,2,3. For a 
system with orthorhombic symmetry in a field along [111] we recover Eq. (j3J), where P is the 
total polarization, a = 2«i, fym] = (4/3)(an+ai 2 ), and qm] = (2/9) (3am +6ai 1 2 +0(123). — 
Experiments in PMN^>22 indicate that X3 > f° r E II [Hl]j implying b\m] < 0. On the 
other hand, for tetragonal symmetry with E || [001] one has 6r oi] = 4«n and c = 6am. For 
example, in PMN— & one finds in this case that % 3 < or 6[ oi] > 0, indicating that «u > 
and a\2 < -an. 



If a(T) > b 2 /c at all temperatures, there is no first order phase transition for E = 0. If 
E > 0, however, a first order phase transition with a jump of polarization AP occurs for E 
exceeding some threshold value E\. This can readily be seen by numerical minimization of 
the free energy (Tjf]) for any pair E, T. As the temperature increases, AP(E, T) decreases and 
vanishes at an isolated critical point T C p, Ecp, 7 '" 22 where the derivative (dP/dT) E diverges.— 
The critical temperature is determined from the equation 



and Ecp from the second of Eqs. ( fT2l) . It should be noted that, generally speaking, b and c 
can also be functions of T.->22 

For E > Eqp the system is in a supercritical regime with continuous temperature and 
field dependence of P(E,T). 

IV. TEMPERATURE AND ELECTRIC FIELD DEPENDENCE OF ECE 

To illustrate the temperature and field dependence of the ECE in relaxors and ferro- 
electrics, we calculate AT(E,T) from Eq. ([7|) for a selected set of parameter values. First, 
we introduce rescaled, dimensionless quantities F and P according to F — > Fv/J and 
P — > P/P max , where P ma x — g/v is the saturation polarization occurring at high field values 
and/or low temperatures. This requires a rescaling of the remaining parameters according 
to a av/(JPl ax ), b bv/(JP^ ax ), c cu/(JI*J, etc. Also, we redefine T -> kT/J 
and E — >• gE/J. From Eq. (fl7|l we see that in relaxors the rescaled parameter a(T) behaves 
as a(T) = T/(l — q) — T , and ai(T) becomes 



Here and until the end of this section, the symbols a, 6, c, ... and P, T, E refer to dimen- 
sionless, rescaled parameters, but elsewhere in this paper the same symbols denote the true, 
physical values of these quantities. 

In ferroelectrics, q = and thus a(T) = T — T , and a± = 1. In relaxors, in the high 
temperature limit T >> 1, q tends asymptotically to zero, thus a ~ T — To and a± ~ 1, i.e., 
the same as in ferroelectrics. 

In Fig. 1, a(T) and ai(T) are plotted for a relaxor ferroelectric with J / J = 0.9 and 
A/J 2 = 0.001. Also shown is the behavior of a normal ferroelectric. It should be noted that 



a C p = a(T C p) 



20~' 



(19) 




(20) 
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the essential difference between ferroelectrics and relaxors is the behavior of the correspond- 
ing coefficients a(T) and a\(T). 

In the following we will choose b = const. = ±1/3 and c = The corresponding ECE 
temperature change is obtained from Eq. (j7j). Using the fact that b\ — c\ — etc., we find 

kT r i 

AT = ^r«i(T) [P(E, T) 2 - P (T) 2 J . (21) 

The polarization P(E, T) will be calculated numerically by simultaneously solving Eq. (jSJ) 
and minimizing the free energy (j3J). We will do that separately for the two cases b > and 
b < 0, assuming the denominator 2vCe to be a constant amplitude factor. 

(i) Case b > 0. As already stated, Pq = P(0,T) = in a relaxor, but in ferroelectrics 
Po(T) 7^ for T < To. The spontaneous polarization Pq(T) is obtained by minimization 
of F(E = 0). In real systems, Pq may not be spatially uniform due to domains. Here we 
assume that \Pq\ has the same value in all domains regardless of their orientation, and that 
the contribution of domain walls to the entropy can be neglected.- 

In Fig. 2(a) we show the calculated values of AT for a relaxor as function of temperature 
for various values of E/Ecp, where Ecp is formally given by Eq. (fl2"]) . although the critical 
point does not exist for b > 0. Also shown in Fig. 2(b) is AT for a ferroelectric with the same 
parameters b, c, but with different a and a\. AT has a peak at T — 0.9 and is in general 
larger than in the relaxor case. At higher temperatures, however, the difference gradually 
disappears. 

(ii) Case b < 0. Eqs. ([HITS]) imply a(0) = ^1 + A/J 2 - J /J = 0.1005. Thus, 
a(0) < 36 2 /16c and the first order phase transition in relaxors at E = is suppressed. 
In a ferroelectric, however, a first order transition in zero field occurs at T\ = 0.9625. The 
critical point is located at Tqp = 1.0398 in relaxors, and at Tqp = 1.015 in ferroelectrics, 
while Ecp = 0.0438 in both cases. In general, AT is found to increase with increasing field 
and exhibits a peak as a function of temperature. In relaxors, the peak position moves to 
higher temperatures with increasing field values, whereas in ferroelectrics the maximum is 
located at T = T 1; where a jump of the spontaneous polarization occurs. 

In Fig. 3, AT(T) is plotted for four values of the field (E/Ecp) 2 , as indicated. At 
E = Ecp in relaxors, there is a jump of AT due to the field- induced first-order transition. 
On the other hand, in the case of a ferroelectric, there is a jump of AT at the zero-field first 
order transition temperature T\ = 0.9625. At the critical point, AT is continuous but with 
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an infinite slope in both cases, while at higher temperatures the difference between relaxors 
and ferroelectrics tends to disappear. 

In Fig. 4, the ECE efficiency AT / E is plotted as a function of E/Eqp for four val- 
ues of temperature close to Tqp- As expected, the maximum efficiency is obtained at the 
corresponding critical points T = T C p and E = E C p- 

Larger values of AT in ferroelectrics rather than in relaxors are mainly due to the sharp 
decrease of cti(T) in relaxors at T < 1 (cf. Fig. 1). However, this does not mean that 
ferroelectrics are better candidates for achieving giant ECE. Namely, one should bear in 
mind that there are other parameters, such as and Ce, which also have a strong impact 
on the ECE. Moreover, the above comparison between relaxors and ferroelectrics makes only 
sense if the coefficients b, c, ... are indeed the same in both cases. Therefore, in discussing 
the ECE in specific systems one should carefully consider the actual physical values of all 
the relevant model parameters, as discussed in the following section. 

V. APPLICATION TO REAL SYSTEMS 

It has been observed experimentally in a variety of systems that the entropy change ASp 
in Eq. (fTUj) is proportional to p 2 } M~M. suggesting that the terms of order ~ P 4 and higher 
in the expansion ([§]) make no contribution. Assuming for simplicity that P = 0, we recover 
from Eq. ([9]) the empirical quadratic relation 



According to the Landau model, in ferroelectrics the partial derivative is equal to 1 or 
= (eo©) -1 at all temperatures, whereas in relaxors it approaches the value ~ 1 at high 
temperatures, but is in general a function of temperature. Thus, in relaxors, /3 is expected 
to be a function of temperature with (3{T) < (e 6) _1 . 

Using the relations ffTUl) and (1221) . the ECE temperature change AT can be written as 



ASp = - X -$P\ 



(22) 



where the coefficient (3 can be expressed through Eqs. (J9j) and ( TT71) . i.e., 




(23) 



AT 



T 



(3P(E,T) 2 . 



(24) 



2C 
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A quadratic relation of the same form is predicted by Eq. (12T|) if the terms bi, ci,... etc. 
are neglected. It is exactly true in the case 61 = c\ — ■ ■ ■ — discussed in Section III. 

Clearly, the above relation is applicable only in the quadratic regime where the empirical 
Eq. (|22p is valid. The parameter /3 can be determined directly from the measured ECE 
temperature change AT, with P(E, T) extracted from dielectric experiments. Alternatively, 
and especially in ferroelectrics or relaxors at temperatures above the freezing temperature, 
we can obtain both P(E, T) and the Curie constant O from the dielectric data, and then 
deduce f3 from Eq. (|23|) . We can then predict AT from Eq. (12^|) . It turns out, however, 
that the values of AT measured directly may differ from the ones deduced from dielectric 
data. For example, in the case of relaxor ferroelectric terpolymer P(VDF-TrFE-CFE) at 
T ~ 300 K and E = 70 MV/M, the ECE measured directly is AT ~ 3.6 K, leading to 
(3 ~ 10 7 V m C _1 K -1 ,— whereas the value deduced from dielectric data is AT ~ 0.87 K.— 
A tentative explanation of this discrepancy is that even far above the freezing temperature 
Tf ~ 277 K the system may still be nonergodic, so that the Maxwell relations, based on 
equilibrium thermodynamics, are not applicable. 30 As a consequence, the value of P(E,T) 
measured on a short time scale is smaller than its thermodynamic long-time limit. Another 
possibility is that the empirical relation Eq. (1221) is an effective quasi-linear relation observed 
in a broad range of large field values, while its derivation based on the Landau expansion is 
by assumption restricted to small fields. 

We can obtain an estimate for the maximum ECE temperature change (AT) max by 
assuming that in a sufficiently strong electric field the polarization reaches its saturation 
value P m ax — g/v. Relation (1221) is not expected to be valid in this saturation regime, and we 
must return to the general expression ([101) . Obviously, the expansion ([9]) cannot be applied 
due to convergence problems. On the other hand, it is well known^ that in the saturation 
regime the excess entropy of the dipolar subsystem tends to zero. Therefore, according to 
Eq. ([8]), Si(P max ) should approach the negative value of the configuration entropy S , i.e., 
Si(P m ax) — > —{k/v) ln(f2). Thus, in the saturation regime Eq. (flOl) leads to 

fcTln(n) 

(Al ) max = — . (25) 

This relation gives the theoretical upper bound on ECE in terms of just three physical 
quantities, v, Ce, and the configuration number Q. Interestingly, this result does not depend 
explicitly on the dipole moment g. Moreover, it does not contain any information about 
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possible phase transition occurring in the quadratic regime. 

For a given value of electric field E, the borderline between the two regimes is expected 
to occur at some temperature T*(E) where the dipolar energy becomes equal to the thermal 
fluctuation energy, i.e., T* = gE/k. For T > T* the system is in the quadratic, and for 
T < T* in the saturation regime. For the above terpolymer P(VDF-TrFE-CFE), we can 
estimate the dipole moment g from the relations P ma x ~ 0.1 C/m 2 and g = k/((3P max ). 
Using the value of /3 determined directly from AT, we have g ~ 5.7 x 10~ 30 C m, and for 
E = 300 MV/m we thus find T* ~ 124 K. 

The average volume associated with a PNR in relaxors, v = V/N, in Eq. (|25p is not 
a priory known and depends on the total number of PNRs. We can estimate v from the 
measured values of and P max using the relation = P max v/ke , and rewrite Eq. (J2"5|) in 
the form 

(AT) ma , - —±-Lpl ax . (26) 

The value of P m ax can be extracted, for example, from hysteresis loops in the saturation 
regime, and from the asymptotic behavior of xi- As already discussed above, this value of 
may differ from the one derived from the experimental value of parameter = (eqQ) -1 , 
observed in ECE experiments in the effective quadratic regime. 

In Table I, the predicted values of (AT) max for a set of selected systems are listed using 
the values of deduced from dielectric experiments. It should be noted that Eqs. (125]) and 
(T26|) provide just a theoretical upper bound for the ECE in the systems listed. In practice, 
the limit of a fully polarized dipolar subsystem might not be accessible because dielectric 
breakdown could occur before complete saturation is reached. Nonetheless, the predicted 
values of (AT) max permit a comparison between various sytems and might be useful in the 
search for a giant ECE. 

Eq. (125]) indicates that a giant ECE is expected to occur in systems with a small value 
of v, or equivalently a large number N of dipolar entities at fixed volume V. For illus- 
tration, let us consider a specific example, i.e., the ferroelectric copolymer P(VDF-TrFE). 
This system consists of microscopic crystalline layers of polarized material embedded in an 
amorphous environment.— Electron irradiation breaks up the layered structure into smaller 
dipolar units, and turns the polymer into a relaxor. Since the number of these new entities is 
now larger than the number of the original microcrystallites, and assuming the same value 
of saturation polarization, one expects a stronger ECE to occur in the irradiated relaxor 
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copolymer than in the original ferroelectric copolymer. This is corroborated by the exper- 
imental values of the coefficient /3 for the irradiated relaxor copolymer P(VDF-TrFE),— 
/3 = 9.3 x 10 7 V m C _1 K _1 , and for the original ferroelectric copolymer in the paraelectric 
phase, 25 = 5.8 x 10 7 V m C^K -1 . 

In the case of inorganic relaxor 8/65/35 PLZT thin films, the value of the coefficient j3 is^ 
~ 1.5 x 10 6 V m C -1 K -1 , whereas for the ferroelectric PZT one finds 1 j3 ~ 7.6 x 10 5 V m 
C -1 K -1 . Again, /3 is larger for the relaxor; however, this comparison seems less conclusive 
since the difference in composition between the two systems is much greater than in the 
above organic case. 

Relation (j25p has been derived here in the framework of Landau theory, however, its 
validity is essentially based on thermodynamic and statistical principles and is hence quite 
general. In particular, it is independent of any mesoscopic models such as the SRBRF model. 
Moreover, it can be easily generalized to magnetic systems, where AT represents the MCE 
temperature change. 

The smallest physical limit for v in ferroelectrics is the volume of the unit cell Vq, yielding 
the ultimate upper bound on (AT) max in Eq. ( 125]) for ferroelectric materials. To emphasize 
this point, let us multiply the numerator and denominator by the Avogadro number iV^ to 
obtain 



Here, R = A^/c is the gas constant and C m = NaVqCe the molar specific heat. This simple 
result does not explicitly contain any information on the microscopic nature of the system, 
i.e., whether it is dielectric, magnetic, etc., since the corresponding electric or magnetic 
dipole moment does not appear in Eq. fl2"Tj) . Therefore, we can regard the above relation 
as a universal law for the theoretical upper bound on the ECE and/or MCE in electrically 
or magnetically polarizable solids. Of course, Eq. ff27|) implies that the polarization or 
magnetization must have reached complete saturation. In the dielectric case, this requires 
large electric fields with the already mentioned possibility of dielectric breakdown. In the 
magnetic case, however, extremely large fields of the order of ~ 100 T may be necessary. 

Eq. f[2"T|) is valid in the saturation regime T < T*, typically at low temperatures, where 
C m is expected to be temperature dependent. In ferroelectrics, the formation of domains 
should be taken into account, and the entropy limit S — > —S could only be reached in 
very high fields. In relaxors, there is an additional difficulty that the relaxation of PNRs at 




(27) 
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low temperatures is extremely slow and the above limit can only be reached after very long 
times. 

At high temperatures, C m normally approaches a certain limit. In metals 
3R according to the Dulong-Petit law, which is generally not valid for complex solids. For 
example, in perovskites, C m ~ 15-R. 10 Thus, the maximum ECE in perovskites at high 
temperatures is of the order (AT/T) 

max r * > 0.0667 ln(f2). This limit applies to a large group of 
perovskites, i.e., paraelectrics, ferroelectrics, relaxors, etc. Specifically, for Q = 8 (see Table 
I), we find (AT/T) 

max ~ 0.139. For example, if the system can support large electric fields 
such that T* > 300 K, this would lead to a giant ECE temperature change (AT) max ~ 40 K 
at room temperature. This exceeds roughly by a factor of ~ 2 the estimated values for 
perovskites in Table I, which were derived from the physical values of the parameters P ma x, 
9, and Ce- 

VI. CONCLUSIONS 

In this paper, we studied the mechanism of ECE in relaxor ferroelectrics (or relaxors) 
and in normal ferroelectrics (referred to as ferroelectrics). Starting from the widely accepted 
result for the ECE temperature change AT(E,T), based on the thermodynamic Maxwell 
relation, we derived an alternative expression which could directly be applied in theoretical 
model calculations. The results for AT(E, T) have been obtained in two physical regimes of 
field-temperature variables E, T. The first of these is the so-called quadratic regime, where 
AT is proportional to the square of the dielectric polarization P(E,T), and the second the 
saturation regime, where P(E,T) is allowed to reach its maximum value Pmax- 

In the quadratic regime, the system can be described by the Landau-type free energy 
model, in which the harmonic Landau coefficient a(T) depends on the physical nature of 
the system, i.e., the behavior of a(T) in relaxors differs from that of ferroelectrics. The 
anharmonic coefficients b, c, which are common to both cases, then determine the critical 
behavior of the system in a field E. The case b > does not show any pronounced anomalies. 
For b < 0, the polarization P and hence AT rises steeply in a relaxor near the isolated critical 
point Ecp, Tcp, whereas in a ferroelectric, the largest effect occurs near the first order phase 
transition in zero field, which is absent in the relaxor case. The ECE efficiency AT jE shows 
a similar behavior in both cases, and at higher fields the same asymptotic values are found. 
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Of course, these conclusions only apply provided that all the remaining physical features, 
such as the number of equilibrium orientations Q of the elementary dipolar entities, etc., are 
the same in both cases. 

Experimentally, the ECE entropy change in the quadratic regime is found to behave as 
AS ~ —(1/2)I3P(E, T) 2 , which is trivially reproduced by the Landau model. The coefficient 
(3 tends to have larger values in relaxors, leading to a stronger ECE. In irradiated organic 
polymer relaxors, this can be explained by the larger number of polar nanoregions (PNRs) 
and thus a smaller average PNR volume v. In the saturation regime, the entropy of the 
dipolar subsystem generally approaches the negative value of the configuration entropy, 
S ~ —(k/v)hi(Q), and the maximum ECE value (AT) max thus crucially depends on the 
orient at ional degeneracy Q. 
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TABLE I: Predicted limiting values of (AT) 

max f° r various perovskite and polymer relaxor systems 

according to Eq. ([26]) . 



Material 


n 


p 

1 max 

(C/m 2 ) 


C E at 300 K 
(10 6 J/m 3 K) 


P = (eoQ)- 1 
(10 5 Vm/C K) 


AT (K) 
at 300 K 


8/65/35 PLZT 


8 


0.65 26 


3.022 


2.4322 


21.4 


PMN-29.5PT 


8 


0.552 


2.8^ 


3.76 7 


24.8 


PMN 


8 


0.55 22 


2.6^ 


3.82M 


20 


P(VDF-TrFE) 68/32 


8 


0.0629 


2.3429 


4.67-29. 


39 



(irradiated film) 
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FIG. 1: Temperature dependence of the coefficient a(T) and its temperature derivative ai(T) for 
a relaxor (solid lines) and a ferroelectric (dashed). 

(a) (b) 




FIG. 2: (a) Calculated temperature dependence of AT in a relaxor with b > and four values of 
electric field E, as indicated, (b) Same, but for a ferroelectric. 
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FIG. 3: (a) Calculated temperature dependence of AT in a relaxor with b < and four values of 
electric field E, plotted on the same vertical scale as in Fig. 1. The critical point lies at Tcp = 1.0398 
and Ecp = 0.0438. (b) Same, but for a ferroelectric with Tqp = 1.05 and Eqp = 0.0438, on 
different vertical scale. Note the first order phase transition at T\ = 0.965. 



(b) 



w 

< 




W 

f- 

< 




FIG. 4: (a) Calculated field dependence of ECE efficiency AT/E in a relaxor with b < 0, at four 
temperatures close to critical temperature Tqp- (b) Same, but for a ferroelectric, with different 
vertical scale. 
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